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We calculate the critical exponents of the chiral phase transition at nonzero tem- 
perature using the thermal and chiral susceptibilities. We show that within a class 
of confining Dyson-Schwinger equation (DSE) models the transition is mean field, 
and that an accurate determination of the critical exponents requires extremely 
small values of the current-quark mass, several order of magnitude smaller than re- 
alistic up- and down-quark masses. In general, rainbow truncation models of QCD 
exhibit mean field exponents as a result of the gap equation's fermion substructure. 



1 Introduction 

It is anticipated that the restoration of chiral symmetry, which accompanies 
the formation of a quark-gluon plasma at nonzero temperature T, is a second- 
order phase transition in QCD with 2 light flavours. Such transitions are 
characterised by two critical exponents: 5), which describe the response of 
the chiral order parameters, X, to changes in T and in the current-quark mass, 
m. Denoting the critical temperature by Tc, and introducing the reduced- 
temperature t := T/Tc — 1 and reduced mass h := m/T, then 

Xoz{-tf , t^Q- , h = 0, (1) 

X (xh^/\ h^0+ , t = 0. (2) 

Calculating the critical exponents is an important goal because of the notion 
of universality, which states that their values depend only on the symmetries 
and dimensions, but not on the microscopic details of the theory. 

The success of the nonlineax. (T-model in describing long-wavelength pion 
dynamics underlies a conjectural that chiral symmetry restoration at finite T 
in 2-flavour QCD is in the same universality class as the 3-dimensional, = 4 
Heisenberg mag net (0(4) model), with critical exponents:El = 0.38 ± .01, 
5^ = 4.82 ± .05. However, recently it was argued that the compositeness 
of QCD's mesons affects the nature of the phase transition and the Gross- 
Neveu model was presented as a counterexample to universality.a Subsequent 
studiesu indicated that nontrivial 1 /N corrections are important and that this 
model has the same critical exponents as the Ising model, as was argued on 
the notion of universality, but only in a scaling, region of width 

Calculating the exponents /3 and S directlyQ from Eqs. (j^) and (||) is often 
difficult because of numerical noise near the critical temperature. Another 
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method is to considerQ the chiral and thermal susceptibihties: 

dX{t,h) 



Xh{t,h) := ^ ' ^ 



dh 



Xtit,h) 



at 



(3) 



At each h, Xi(i, h), i = h,t, are smooth functions of t with maxima Xi'^ the 
pseudocritical points tp^. Near the critical point t — — h we have 

h^/^'^'Wtl,, (4) 

xT = Xh{t^c,h) cx 2„:=l-i, (5) 

xr = Xt(4c, « h-^' , ^ (1 - /3) . (6) 

Therefore, by calculating the chiral and thermal susceptibilities and locati|ii|g 
the pseudocritical points, one can determine Tc and the critical exponents.EliD 

2 Quark Dyson Schwinger Equation 

We have analysedi Xh{t,h) and Xt(i, /i) in a class of confining DSE models 
that underlies rriany successful phenomenological applications^ at both zero 
and finite- (T, The foundation of our study is the renormalised quark DSE 

S^'^{p,uJk) ■■= ■ pA{p'^,u!k) + i-iiUJk C[p^,ujk) + B[p^,ujk) (7) 

Here LUk — (2fc + 1) ttT is the fermion Matsubara frequency and mR{C,) is the 
current quark mass at the renormalisation point C,. The self-energy is 

S'(p,t^fe)=r^ j j—^ -g'^ D^^{p- q,^k-i)lt,S{q,uJi)T^ , (9) 

/ — — oo 

with D^^{k, flj) the renormalised dressed-gluon propagator and the renor- 
malised dressed-quark-gluon vertex. In renormalising the DSE we require 

•S^^Ip, wo)|p2+^2^^2 ■p + ij4(^o+mR . (10) 

Equations (^-(|l0|) define the exact QCD gap equation. 

We use the rainbow trunctation for the vertex, Ti, — "fi^, which is the lead- 
ing term in a 1/iVc-expansion of the vertex, and consider three models in which 
the long-range part of the interaction is an integrable infrared singularity,tJ 
motivated by T = studies of the gluon DSE:llil 

g'D^^ik, n,) = P^,{k, n,)V{k, n,;mg) + pj;(fc, n,)V{k, n,;0) , 
'D{k,rtj;mg) ■.= 2TT^D^5ojS^{k)+VM{k^ + n^j+ml), (11) 
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where Pl^ = Pi = 0,P^,= S,, - hkj/k\ P^, = S,,, ~ k^K/ik^ + n^) - Pj,, 
iTig is a Debye mass, and Z? is a mass-parameter fitted to m^r and /jr at T = 0. 
We compare the results for 3 different models, denoted by Vm, M = A, B, C. 

One order parameter for dynamical chiral symmetry breaking is the quark 
condensatet3 (99)^- There are other, equivalent order parameters and in cal- 
culating the chiral and thermal susceptibilities we employ 



^0,^0), ^c:=^^^^^y (12) 



They should be equivalent and, as we will see, the onset of that equivalence 
is a good way to determine the /i-domain on which Eqs. (^)-(^) are valid. 
Further, we have verified numerically that in the chiral limit (m = 0) and for 
t ^ 0: /tt oc (qq) cx X{t,0); i.e., that these quantities are all equivalent, bona 
fide order parameters. It thus follows from the pseudoscalarjjuass formulail3 
f^m^ — 2mii{(^){qq)^, that increases with temperature.El 

3 Results 

The first model we consider is an infrared, dominant model with I?a(s) = 0, 
and the mass-scale D — 0.56 GeV^ fixedE2l by fitting tt- and p-meson masses 
at T = 0. A current-quark mass of m = 12 MeV yields = 140 MeV. The 
quark DSE obtained with Va is an algebraic equation. Chiral symmetry 
restoration is therefore easy to analyse and either directly, via Eqs. (|l|) and 
(^, or using the susceptibilities and Eqs. it is straightforward to 

establish^ that this model has mean field critical exponents and to determine 
the critical temperature in Table 0. The exponents are unchangecO and Tc 
reduced by < 2% upon the inclipion of some higher-order 1 /A^c-corrections to 
the dressed-quark-gluon vertex.E3 

3.1 Model B: QED-like tail 

To improve the ultraviolet behaviour, we consider a model! with 

ifi o 1 - e~'*/("^'"*) 
2?b(,) = H,2^^ , (13) 

9 S 

and D = (8/9) m^. The mass-scale rrit = 0.69 GeV = 1/0.29 fm marks the 
boundary between the perturbative and nonperturbative domains, and was 
fixedtil by requiring a good description of tt- and p-meson properties at T = 0. 

The quark DSE obtained with this model can be solved numerically and 
Xh^{h) and Xt"^(^) ^re depicted in Fig. |l|(a). Following Eqs. (||) and (||), the 
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h iog,„(h) 

Figure 1. (a) yS^{h) and x^'^i^) calculated in model B. The slope of the straight lines is 
given m Table m and tHey are drawn through the smallest h- values; (b) Zf^{h) (top) and zt{h) 
(bottom) from Eq. (|l4[). The dashed lines are the mean field values: Zf^ = 2/3, zt = 1/3. 



critical exponents can be determined by defining a local critical exponent as 
a function of the reduced mass h for each of the equivalent order parameters: 

z.{h):^ ~h^^, (14) 

see Fig. ^b). h lies in the scaling region when Zi is independent of the order 
parameter. This shows that the scaling relations are not valid until 

logio(V^n) < -7, (15) 

where hu = nif/jTc corresponds to the current-quark mass that gives 
TO^ = 140 MeV in this model. The values of Zh and zt in Table |l| are obtained 
by a Pade fit to the five smallest /i- values in Fig. 0(b), and extrapolating to 
/i — !■ 0+. The critical temperature is obtained using Eq. (Q); its value is in- 
sensitive to whether t^^ or t^^ is used and to which of the equivalent order 
parameters is used. 





mean field 


A 


B 


C 


Te (MeV) 




169 


174 


120 


Zh 


2/3 


0.666 


0.67 ± 0.01 


0.667 ± 0.001 


Zt 


1/3 


0.335 


0.33 ± 0.02 


0.333 ± 0.001 



Table 1. Critical temperature for chiral symmetry restoration and critical exponents char- 
acterising the second-order transition in the three exemplary models. 
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3.2 Model C: Logarithmic tail 

Finally, we consider the finite-T extension§ of a model which further improves 
the ultraviolet behaviour, via the inclusion of the one-loop In-suppression at 
s ^ Aqcd- Again, the parameters are fixed at T = by requiring a good 
fit to a range of tt-, K-meson properties;0 recent xialculations show that the 
vector mesons are also described well in this modelJij The study of chiral sym- 
metry restoration in this model is very similar to the previous study, with the 
additional In [s] -suppression in the ultraviolet making the numerical analysis 
easier. The critical temperature and exponents are presented in Table ^ Also 
in this case the scaling relations are only valid for very small current-quark 
masses: \ogif^{h/hu) < —5. These results are qualitatively, and for the crit- 
ical exponents quantitatively, independent of the parameters in this model. 
Direct calculation of the critical exponents using Eqs. (|l|) and (Q) are in good 
agreement with the critical exponents found using the susceptibilities. 



4 Conclusions 

It is clear from Table |l| that each of these models is mean field in nature. In 
hindsight that may be not surprising because the long-range part of the in- 
teraction is identical. However, the models differ by the manner in which the 
interactions approach their long-range limits, and our numerical demonstra- 
tion of their equivalence required extremely small values of the current-quark 
mass, Eq. (^5|). This might also be true in QCD; i.e., while Tc is relatively 
easy to determine, very small current-quark masses may be necessary to accu- 
rately calculate the critical exponents from the susceptibilities. In that case, 
calculation of (3 and 6 via lattice-OpD will not be easy. The discrepancies 
found in recent lattice calculationsta could be a signal of this difficulty. 

The class of models we haye considered can describe the long-wavelength 
dynamics of QCD very wellBu in terms of mesons that are quark-antiquark 
composites. The characteristic feature is the behaviour of the confining in- 
teraction. It provides a driving term in the quark DSE proportional to the 
dressed-quark propagator, which means that boson Matsubara zero-modes do 
not influence the critical behappur determined from the gap equation. The 
class of Coulomb gauge modelfO also describes mesons as composite particles 
and it too exhibits mean field critical exponents. The long-range part of the 
interaction in that class of models corresponds to the regularised Fourier am- 
plitude of a linearly rising potential. Hence it is not equivalent to ours in any 
simple way, except insofar as zero modes do not influence the gap equation. 

The quark DSE is the QCD gap equation and the many equivalent chiral 
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order parameters are directly related to properties of its solution. We have 
observed that several classes of models exhibit the same (mean field) critical 
exponents. Only in our simplest confining model did we consider the effect 
of f /A^c-corrections to the quark-gluon vertex, and in that case the critical 
exponents were unchanged. These results suggest that mean field exponents 
are a feature of the essential fermion substructure in the gap equation. It 
can likely only be false if nonperturbative corrections to the vertex are large 
in the vicinity of the transition. In this context the role of mesonic bound 
states, which can appear as nonperturbative contributions in the dressed- 
quark-gluon vertex, has to be studied in more detail. This might also give a 
nontrivial dependence on the number of fermion flavours, as is anticipated on 
universality arguments. 
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